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We present a symmetry-based theory of the depolarization induced by subwavelength metal hole 
arrays. We derive the Mueller matrices of hole arrays with various symmetries (in particular square 
and hexagonal) when illuminated by a finite-diameter (e.g. gaussian) beam. The depolarization is 
due to a combination of two factors: (i) propagation of surface plasmons along the surface of the 
array, (ii) a spread of wave vectors in the incident beam. 

PACS numbers: 78.67.-n, 73.20.Mf, 42.25.Ja, 02.20.-a 



The transmission of a metal film perforated with a two- 
dimensional periodic array of subwavelength holes is ex- 
traordinary large due to the resonant excitation of surface 
plasmons (SPs) ^ . Recently, the polarization properties 
of this enhanced transmission have drawn considerable 
attention 0jS0i0,0i0|- ^ particularly intriguing polar- 
ization property is that illuminating the array by a fully 
polarized beam at normal incidence leads to important 
reduction of the degree of polarization in transmission, 
corresponding to a strongly space- variant character of the 
polarization state of the output beam . This seems in 
clear contrast with what could be guessed from general 
symmetry arguments. For plane wave illumination, a pe- 
riodic structure can theoretically only transform a pure 
input state of polarization (SOP) to another pure out- 
put SOP (see below). Thus, the observed depolarization 
must be critically related to additional spatial degrees of 
freedom of the system. This is exactly what excitation of 
SPs can potentially yield since SPs propagate in distinct 
directions related to the symmetry of the hole array. 

The aim of this Letter is to provide a proper theo- 
retical framework describing these depolarization effects. 
Formally, spatial symmetries of arrays within the context 
of SP excitations must be confronted with photon polar- 
ization SU(2) symmetry, a context where Mueller algebra 
is the natural tool |^ ■ The central issue of our work 
is that implementations of symmetry operations have to 
account for the fact that SPs are a source of effective 
spatial dispersion at the arrays interfaces but only so if 
they are excited by a finite-diameter beam (e.g. a gaus- 
sian beam). As we will discuss, it is the strong angular 
dependence related to such a spatially dispersive trans- 
mission that can lead to sizeable depolarization. 

Any two-dimensional pattern can only match five pos- 
sible lattices: the five 2D Bravais lattices 10] . In order 
to complete such patterns into 2D crystal structures, a 
primitive cell is associated with each lattice point. The 
required symmetry compatibility of the primitive cell and 
the Bravais lattice leads to the well-known seventeen 2D 
point groups. Here we restrict ourselves to circular holes 
on a Bravais lattice, in other words to the simplest prim- 
itive cell with full symmetry. Therefore, the point group 



symmetry of our arrays is reduced to the spatial symme- 
try of the chosen Bravais lattice which is a Cnv group, 
where n denotes rotations by 27r/n radians about the 
origin (with n = 1,2,3,4 or 6) and v refers to mirror 
symmetries (if allowed). 

The optical transmission of a given array under gen- 
eral conditions of illumination is fully described by the 
transmission matrix which relates, at a specific wave- 
length A, the transmitted field E°"' to the incoming field 
E'". In the paraxial approximation used throughout this 
Letter, far field angles {9x,9y), collected in a 2-column 
vector 9, define paraxial fields as two independent com- 
plex numbers forming a spinor E = (i?i, E2) in the cho- 
sen (l,2)-basis. The transmission matrix t{X;9) of the 
zero-order diffracted beam in the far field is therefore a 
2x2 matrix 



E°"' (A; 6*) =t(A;e') E'"(A;6') 



(1) 



From this input-output description, spatial symmetries of 
the chosen array define orthogonal transformation ma- 
trices B that leave the transmission matrix unchanged 
via BH (A; 9) B — t (A; B9). For the case of plane wave 
illumination at arbitrary angle 9, spatial symmetry op- 
erations therefore merely relate transmission matrix ele- 
ments with various angular arguments. It is only when 
one resorts to angularly integrated expressions, i.e. to 
illumination by a rotationally symmetric beam (as for 
instance a gaussian beam at normal incidence), that one 
recovers simple symmetry-based relations between input 
and output SOPs of a given beam. 

Depolarization effects are generally most efficiently 
addressed when one resorts to the Stokes parameters 
(So, Si, S2, S3). These parameters, together with the use 
of the associated Mueller algebra, define the natural the- 
oretical habitat of polarization properties in connection 
with symmetry arguments 01. The Stokes parame- 
ters are real-valued and represent four intensity measure- 
ments on a light beam: So = (/) corresponds to the total 
intensity in the beam, as measured without any polariza- 
tion selection, and Si = (/o° — Igo°), S2 = (/450 — 
S3 = — I„- ) represent three balanced intensity mea- 
surements where 0°, 90° and ±45° refer to orientations of 
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FIG. 1: Cartoon-like representation of a hole array, chosen 
here as a square array with circular holes, illuminated by 
a finite-diameter (gaussian) beam. The input beam, chosen 
here as a pure horizontal (0°) SOP, excites SPs along the 
principal diagonals of the array (note that the reciprocal lat- 
tice coincides with the direct lattice). On the backside of 
the array, the corresponding near field pattern is shown with 
SPs propagating and polarized along the -1-45°/ — 135° and 
-1-135°/ — 45° lobes, onto which the input SOP can thus be 
projected. The central region corresponding to the incident 
beam remains polarization isotropic but when moving radially 
outward, the input horizontal polarization is modified through 
elliptically polarized states towards ±45°-linear polarizations 
due to the longitudinal character of the SPs. As symbol- 
ized in the figure, the output SOP becomes space- variant, the 
variation in polarization states being parametrized by the far 
field angles 9 — {9x,0y)- This 2D far field pattern is Fourier 
related to the sketched near field pattern with polarization 
directions remaining unmodified. Therefore, the polarization 
directions observed in the far field, and more particularly at 
large angles, are (maybe somewhat surprisingly) associated 
with directions along which corresponding SPs do not propa- 
gate, since a wide near field corresponds to a narrow far field. 
The same general picture remains valid with an input circular 
SOP. In this case, the input SOP passes through various el- 
liptically polarized states to reach the linear SOPs along the 
diagonals. Extension of this description to an hexagonal array 
is straightforward. 



a linear analyzer and cr^((T^) to a right (left)-handed cir- 
cular analyzer. The brackets (• ■ •) stand for averaging 
over the spread of the wave vectors of the beam, i.e. an 
integration over the far field angles (see Fig. 1). A simi- 
lar integration in time has no effect, as the range of input 
frequencies is assumed to be much smaller than the SP 
resonance structure. The degree of polarization 11 of the 
beam is given by 11 = [{Sj + + S§)/Sg]i/^. For unpo- 
larized light, 81,82,83 components are "averaged-out" 
by the angular integration so that 11 = 0. For a fully 
polarized beam, 11 = 1 whereas for partially polarized 
light, < n < 1. 

Collecting the 8tokes parameters in a four-vector S = 
(80,81,82,83), transformation by the array from input 
to output 8tokes vectors S°"* = MS'" is given by the 



4x4 real-valued Mueller matrix M when assuming an 
arbitrary input 80P with 11 = 1. This matrix completely 
describes all possible changes of this initially pure 80P 
due to transmission through the array; it is directly re- 
lated to the transmission matrix t appearing in Eq. (1). 
The global structure of the M matrix can be derived after 
we expand input and output Stokes vectors in terms of 
SU(2) Pauli matrices generators (ctq, <^i, "'2, fa). This ex- 
pansion fixes the SU(2) decomposition of Mueller matrix 
components to 

Mq/3 = - {Vijkl{Ca)ji{(^f3)kl) , (2) 

where we used summation over repeated indices. This 
decomposition is based on the components Vijki — 

Uk (A; e) t*i (A; 6*)^ of the 4'h-rank tensor T' = (t ® t*), 
the star symbol denoting complex conjugation and the 
"circled-cross" symbol a tensorial product on t-matrices. 
Symmetry properties of the Mueller matrix will result 
from conjunction between spatial symmetries of the lat- 
tice, as implemented at the level of the averaged V com- 
ponents, and the 8U(2) symmetry of Pauli matrices. We 
assume a fully polarized input beam, i.e. 11'" = 1 (see 
below for a discussion of the general case), so that the 
average is performed over the far field angles present in 
the input beam. More precisely, this average corresponds 
to an intensity- weighted integration over the input angles 
via {■■■) = J d9---r''{X;9). 

The simple case of plane wave illumination of the ar- 
ray at normal incidence only retains the 6 = value in 
the averaging process. It is then straightforward to show 
that for square {C4v) and hexagonal (Cev) spatial symme- 
tries, both the transmission matrix and, as immediately 
seen from Eq. (2) , the Mueller matrix are proportional to 
the identity matrix. So we find, not surprisingly, perfect 
polarization preservation for transmission through such 
arrays. For rectangular {C2v) symmetry, the SOP may 
change (corresponding to birefringence and/or dichro- 
ism) but the degree of polarization is conserved. The 
situation becomes fundamentally different when the ar- 
ray is illuminated with a finite-diameter beam and when, 
at the same time, SPs are excited at the array interfaces. 
Eq. (2) implies that these two aspects have to be com- 
bined in the transmission process, as we explain now. 

First, the finite-diameter of the incident beam corre- 
sponds to a spread of incoming wave vectors which calls 
for the averaging in the definition of the Stokes param- 
eters and the V tensor. Second, SP excitations stick to 
preferred propagation directions on the array interfaces 
according to the well-known dispersion relation which re- 
stricts the propagation of SP waves on the reciprocal lat- 
tice of the array 'l l'- SP excitations are then possibly de- 
localized beyond the excitation spot of the incident beam, 
the condition being that the transverse coherence length 
of the input beam at the surface of the array should be 
smaller than the SP coherence length, which can among 
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others be estimated from spectral widths of transmission 
peaks associated to SP resonances Q . Identifiable propa- 
gation axes for SPs induce an effective spatially dispersive 
response of the illuminated array; this is analogous to the 
occurence of spatial dispersion in crystal optics |0| , the 
hole array acting as a 2D crystal. 

In this situation, a space- variant output SOP is read- 
ily expected. Since the SPs are dominantly longitudinal 
surface waves, the excited areas on the hole array ex- 
tending beyond the excitation spot are mainly linearly 
polarized along the axes of the reciprocal lattice assigned 
to the excited SP mode. Polarization isotropy is found 
only in the excitation spot itself whereas a progressive 
modification of the input polarization is involved in the 
propagation of SPs, towards linear polarizations oriented 
along axes of the reciprocal lattice. In the Fourier-related 
far field, as sketched in Fig. (1), the output SOP con- 
tains therefore various angularly separated polarization 
components, i.e. it is (partly) depolarized. Most essen- 
tial for our analysis is that, due to the far field angle 
dependence of the transmission matrix amplitudes tij, 
symmetry operations of the chosen array must be car- 
ried out within the angular integration, bringing forth 
different constraints on V for the case of a rotationally 
symmetric beam versus a plane wave as input illumina- 
tion. 

More specifically, spatial symmetry operations are 
easily analyzed by referring to the circular-state basis 
(cr+,cr_). Mirror symmetry exchanges tT+ into cr_ and 
vice- versa and is directly related to complex conjugation. 
Rotation simply corresponds to phase factors induced 
on transformed elements. As a 4"^-rank tensor, V fea- 
tures in this basis a priori 16 independent components. 
Mirror symmetry in combination with angular integra- 
tion reduces this to 8 independent ones and moreover 
implies that these 8 components are real-valued. Then, 
four- and six-fold rotational invariance retain only 4 in- 
dependent components for the former and 3 for the lat- 
ter. Restricting the V tensor to these 4 and 3 compo- 
nents, rotational symmetries together with mirror sym- 
metry force the ctq and CT/j components to the tight Sap 
"selection-rule". This entails that Mueller matrices for 
square and hexagonal arrays must be diagonal, with the 4 
on-diagonal elements given in the (a;, y) laboratory frame 
as 

Mo = {\txx\^ + \tyy\^ + + \tyx\^) 

Ml = {\t..\^ + \tyy\^~\t.y\^~\ty.\^) 
M2 = 2n[{t^jly+t,ytl^)] 

M3 = 25R [{t^Jly ~ t^yt;,)] , (3) 

omitting A and 9 dependences. In practice, the trans- 
mission matrix components tij depend on the actual 
structure of the chosen array, in particular on the 
dimensions of the holes and on the permittivity of the 
metal used (real and imaginary parts). These param- 



eters characterize in fact the SP modes on the array 
interfaces; their quantitative evaluation (i.e. beyond 
symmetry aspects) requires extensive model calculations 
that are outside the scope of this Letter. As special 
cases, we find that the Mueller matrix M is proportional 
to the identity matrix both in the limit of zero hole size 
and that of highly dissipative metal permittivity (for 
e.g. Cr) where SPs cannot be excited. 

The diagonal forms of the Mueller matrices for square 
and hexagonal arrays immediately reveal that there is 
no mixing of Stokes parameters so that the excitation 
does not produce any new spatially averaged polariza- 
tion in the transmitted beam. With no enpolarizing 
capabilities ;13i|. these arrays act as purely depolarizing 
optical elements. To each input Stokes parameter 
i = (1,3) is associated a given degree of polarization 
-Qout _ Mi/Mo < 1, which depends on the opening 
angle of the input beam. For an hexagonal array, the 
reduction from 4 to 3 independent V components results 
in 3 independent elements of the Mueller matrix with 
Ml — M2. As a consequence, two input beams with 
the different pure SOPs S'"/S|," = (1, ±1,0,0) and 
Si"/Sj," = (1,0, ±1,0) wiU suffer the same amount of 
depolarization in transmission. Due to the absence 
of mixing of Stokes parameters, this can be measured 
directly through standard crossed polarization analysis. 
It is easy to show that such a direct measurement of 
depolarization is no longer possible for a rectangular 
array due to the loss of the simple diagonal structure 
of the Mueller matrix when going to two-fold symmetry 
0. 

In conclusion, we have reported a symmetry-based 
theory of the optical polarization properties of nanofab- 
ricated metal hole arrays in the presence of SPs excited 
on the arrays interfaces by a finite-diameter beam. Our 
theory has two essential ingredients: (i) multi-plane 
wave illumination and (ii) spatial dispersion related 
to SP excitations. The former corresponds to a far 
field angular integration performed at the level of V 
components, and the latter implies that this integration 
is done over transmission amplitudes depending actually 
on far field angles. A tensorial symmetry analysis shows 
that such arrays induce depolarization in transmission, 
a fundamental property which has not been addressed 
before. We stress that our analysis is based on the 
assumption of a fully polarized input SOP. If the input 
beam itself is already partially depolarized, that is if 
11'" < 1, the angular integration cannot be performed 
separately on the Mueller matrix M and on the input 
Stokes vector S"', bringing additional complexity into 
the problem. Additionally, and in view of the recent 
discovery of SP-assisted polarization entanglement 2] , it 
would be very interesting to develop a quantum version 
of our approach which could be based upon twin-photon 
Stokes parameters 15]. 

As a further outlook, we have started a theoretical 



and experimental extension of the present work to 
depolarization effects with non-circular holes. Also 
quasipcriodic structures, such as 2D Penrose quasicrys- 
tals [13, are of particular interest; there, the appearance 
of forbidden point group symmetries should lead to 
unexpected depolarization behaviours. 
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